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Abstract 

We calculate transverse momentum distribution of direct photons from various sources by taking into account the 
initial state momentum anisotropy of quark gluon plasma (QGP). The total photon yield is then compared with the re- 
cent measurement of photon transverse momentum distribution by the PHENIX collaboration. It is also demonstrated 
that the presence of such an anisotropy can describe the PHENIX photon data better than the isotropic case in the 
present model. We show that the isotropization time thus extracted lies within the range 1 .5 > t 1so > 0.5 fm/c for the 
initial condition used here. 
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1. Introduction 

The primary goal of relativistic heavy ion collisions is to create a new state of matter, called quark gluon plasma 
and to study its properties through various indirect probes. Out of all the properties of the QGP, the most difficult 
problem lies in the determination of isotropization and thermalization time scales (t[ so and r t h e i- m )- Studies on elliptic 
flow (upto about pj ~ 1 .5 - 2 GeV) using ideal hydrodynamics indicate that the matter produced in such collisions 
becomes isotropic with t; so ~ 0.6 fm/c ( 1 ). On the other hand, using second order transport coefficients with conformal 
symmetry it is found that the isotropization/thermalization time has sizable uncertainties (2). Consequently, there are 
uncertainties in the initial temperature as well. Electromagnetic probes have been proposed to be one of the most 
promising tools to characterize the initial state of the collisions (3, 4). Because of the very nature of their interactions 
with the constituents of the system they tend to leave the system without much change of their energy and momentum. 
In fact, photons (dilepton as well) can be used to determine the initial temperature, or equivalently the equilibration 
time. 

It is to be noted that while estimating photons from QGP (|5j|6lE]l, it is assumed that the matter formed in the rela- 
tivistic heavy ion collisions is in thermal equilibrium. The measurement of elliptic flow parameter and its theoretical 
explanation also support this assumption. On the contrary, perturbative estimation suggests the slower thermalization 
of QGP ©. However, recent hydrodynamical studies (2) have shown that due to the poor knowledge of the initial con- 
ditions there is a sizable amount of uncertainty in the estimate of thermalization or isotropization time. In view of the 
absence of a theoretical proof behind the rapid thermalization and the uncertainties in the hydrodynamical fits of ex- 
perimental data, such an assumption may not be justified. Hence in stead of equating the thermalization/isotropization 
time to the QGP formation time, in this work, we will introduce an intermediate time scale (isotropization time, Tj so ) to 
study the effects of early time momentum-space anisotropy on the total photon yield and compare it with the PHENIX 
photon data d9l fTQl fTTb . Recently, it has been shown in Ref. (fl2l) that for fixed initial conditions, the introduction of 
a pre-equilibrium momentum-space anisotropy enhances high energy dileptons by an order of magnitude. In case of 
photon transverse momentum distribution similar results have been reported for various evolution scenarios (13). 
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The plan of the paper is the following. In the next section we will discuss the mechanisms of photon production 
from various possible sources and the space-time evolution of the matter very briefly. Section 3 is devoted to describe 
the results for various initial conditions and we summarize in section 4. 

2. Formalism 

2.1. Photon rate: Anisotropic QGP 

The lowest order processes for photon emission from QGP are the Compton (q(q) g — > q(q) y) and the annihilation 
(qq — > gy) processes. The rate of photon production from anisotropic plasma due to these processes has been 
calculated in Ref. (TT4l . The soft contribution is calculated by evaluating the photon polarization tensor for an oblate 
momentum-space anisotropy of the system where the cut-off scale is fixed at k c ~ y[gpi W rd- Here phard is a hard- 
momentum scale that appears in the distribution functions. The differential photon production rate for 1 + 2 — > 3 + y 
processes in an anisotropic medium is given by (fT4l) : 

p dN _ N r d 3 Pl d 3 p 2 d 3 p 3 
d 4 xd 3 p ~ 2(2tt)3 J 2£!(2/r)3 2E 2 (2tt) 3 2E i (2nf ™' Pbari ' & 

X (27T) 4 S( Pl +P2- P 3 ~ P)\M\ 2 [1 ± / 3 (p 3 , .Phard, £)] (D 

where, \Ai\ 2 represents the spin averaged matrix element squared for one of those processes which contributes to the 
photon rate and N is the degeneracy factor of the corresponding process. £ is a parameter controlling the strength of 
the anisotropy with £ > -1. fufi and fj are the anisotropic distribution functions of the medium partons. Here it is 
assumed that the infrared singularities can be shielded by the thermal masses for the participating partons. This is a 
good approximation at short times compared to the time scale when plasma instabilities start to play an important role. 
The anisotropic distribution function can be obtained ( 15 ) by squeezing or stretching an arbitrary isotropic distribution 
function along the preferred direction in momentum space, yj-(k, Phard) = f' w ( V* 2 + £(k.n) 2 , pi W rd), where n is the 
direction of anisotropy. It is important to notice that £ > corresponds to a contraction of the distribution function in 
the direction of anisotropy and - 1 < £ < corresponds to a stretching in the direction of anisotropy. In the context 
of relativistic heavy ion collisions, one can identify the direction of anisotropy with the beam axis along which the 
system expands initially. The hard momentum scale phard is directly related to the average momentum of the partons. 
In the case of an isotropic QGP, pi, ar d can be identified with the plasma temperature (T). 

2.2. Photon rate : Isotropic case 

As mentioned earlier the QGP evolves hydrodynamically from t; so onwards. In such case the distribution func- 
tions become Fermi-Dirac or Bose-Einstein distributions. The photon emission rate, in isotropic case, from Compton 
(<7@ g ~ * "7® j) an d annihilation (qq — > gy) processes has been calculated from the imaginary part of the photon 
self-energy by Kapusta et al. (fT6t in the 1-loop approximation. However, it has been shown by Auranche et al. dTTT l 
that the two loop contribution is of the same order as the one loop due to the shielding of infra-red singularities. 
The complete calculation upto two loop was done by Arnold et al. (18). In this paper we have calculated the photon 
production rate from hot hadronic matter. We follow the calculations done in Ref. ( 19 ) where convenient parameteri- 
zations have been given for the reactions considered. These parameterizations will be used while doing the space-time 
evolution to calculate the photon yield from meson-meson reactions. The photon emission rate (static) from reactions 
of the type B M — » By (B denotes baryon) has been calculated in Ref. d20b . It is shown that this contribution is 
not negligible compared to that meson-meson reactions. To evaluate photon rate due to nucleon (and antinucleon) 
scattering from n, p, cj, tj and a\ mesons in the thermal bath we use the phenomenological interactions described in 
Ref. (20). Besides the thermal photons from QGP and hadronic matter we also calculate photons from initial hard 
scattering from the reaction of the type Ha hs — > yX using perturbative QCD. We include the transverse momentum 
broadening in the initial state partons (I2TII221 . 
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Figure 1 : (Color online) Medium photon spectrum, dN/d 2 prdy, at y = for the free-streaming interpolating model (6 = 2) for three different 
values of isotropization time, r, JO , with initial conditions, T\ = 440 Gev and Ti = 0.1 fm/c. 



2.3. Space-time evolution 

The expected total photon rate must be convoluted with the space-time evolution of the fireball. The system 
evolves anisotropically from r, to T[ so where one needs to know the time dependence of p^-d and We have used a 
phenomenological model (fl2l PPSl ) to describe the time dependence of phard an d f ■ In the frame work of this model, 
£ = at t = Ti and it grows with time (t) and reaches maximum at t = t,- jo , after that £ decreases to zero at r >> r, so . 
We shall follow the work of Ref. JT2l [131 to evaluate the p T distribution of photons from the first few Fermi of the 
plasma evolution. In our calculation, we assume a first-order phase transition beginning at the time T c iPhmi(jc) = T c ) 
and ending at r# = rdT c where rj = gQlgH is the ratio of the degrees of freedom in the two (QGP phase and hadronic 
phase) phases. Therefore, the total thermal photon yield, arising from the present scenario is given by, 



dN 
cP-pjdy 



I 



d 4 xE 



dR 

~¥p 



I 



d xE 



dR 

d^p 



(2) 



hydro 



where the first term denotes the contribution from the anisotropic QGP phase and the second term represents the 
contributions evaluated in ideal hydrodynamics scenario. 



3. Results 

We have considered the initial condition, Ti = 440 MeV, t,- =0.1 fm/c and free-streaming interpolating model 
(5 = 2) (TPJl l23T l for the pre-equilibrium evolution. In this initial condition the maximum value of £ will be ~ 70 
at t = t, so . In Fig. ([T]l we present the photon yield due to Compton and annihilation processes in the mid rapidity 
(8 y = 7r/2, y being the angle between the photon momentum and the anisotropy direction) as a function of photon 
transverse momentum. In estimating this result, we have used a s = 0.3. Different lines in Fig. [T] correspond to 
different isotropization times, t; so . We clearly observe enhancement of photon yield when r; so > t, . The enhancement 
of photon yield in the transverse directions (y = 0) is due to the fact that momentum-space anisotropy enhances the 
density of plasma partons moving at the mid rapidity ( fT3l) . To show that the presence of initial state momentum 
anisotropy and the importance of the contribution from baryon-meson reactions we plot the the total photon yield 
assuming hydrodynamic evolution from the very begining as well as with finite r; so (right panel describes the total 
contribution with and without the initial state momentum space anisotropy only for Tj so = 1 fm/c) in Fig. (j2j). It is 
clearly seen that some amount of anisotropy is needed to reproduce the data. We note that the value of t; so needed to 
describe the data also lies in the range 1.5 fm/c> Tj so > 0.5 fm/c for both values of the transition temperatures. 



Lusaka Bhattacharya et al. /Nuclear Physics A 00 (2011 ) 1^ 



4 



i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — r 




p T (GeV/c) p T (GeV/c) 



Figure 2: (Color online) Photon p T distributions at RHIC energies with initial condition T, = 440 GeV, for (a) T c = 192 MeV and (b) 170 MeV. 



4. Conclusion 

To summarize, we have calculated total single photon transverse momentum distributions by taking into account 
the effects of the pre-equilibrium momentum space anisotropy of the QGP and late stage transverse expansion on 
photons from hadronic matter with various initial conditions. To describe space-time evolution in the very early stage 
we have used the phenomenological model described in Ref. (12) for the time dependence of the hard momentum 
scale {phard) and plasma anisotropy parameter (£). To calculate the hard photon contributions we include the transverse 
momentum broadening in the initial hard scattering. The total photon yield is then compared with the PHENIX photon 
data. Within the ambit of the present model it is shown that the data can be described quite well if Ti s0 is in the range 
of 0.5 - 1.5 fm/c for all the combinations of initial conditions and transition temperatures considered here. It is to 
be noted that the apparent hump observed in all the figures (except Fig. (5)) needs to be understood and we wish to 
discuss it in a subsequent paper. 
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